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Abstract 



. In this paper we study the work of Leinweber by applying the Continuum Model of 

I QCD Sum Rules (QCDSR) to the analysis of (quenched) lattice correlation functions. 

We expand upon his work in several areas: we study meson states as well as baryons; we 
analyse data from several lattice spacings; and we include data from the Sheikholeslami- 
Qh! Wohlert (clover) improved action. We find that the QCDSR Continuum Model Ansatz 

can reproduce the data, but only for non-physical values of its parameters. This leads 
us to reject it as a model for hadronic correlation functions. 

We study the non-relativistic quark model and conclude that it predicts essentially 
the same form for the correlation function as the QCDSR Continuum Model approach. 
5^ I Furthermore, because it doesn't have the Continuum Model's restrictions on the pa- 

rameters, the non-relativistic quark model can be viewed as a successful Ansatz. 

As well as studying the validity or otherwise of the QCDSR Continuum Model ap- 
proach, this paper defines 4-parameter fitting functions that can be used to fit lattice 
data even for a time window close to the source. These functions are shown to be an 
improvement over 2-exponential fits especially in the case of mesons. We encourage the 
application of this approach to situations where the conventional fitting procedures are 
problematic due to poor ground state dominance. 
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1 Introduction 



Lattice Gauge Tlieory and QCD Sum Rules (QCDSR) are two areas of researcli 
wliich liave been widely employed to deepen the understanding of systems gov- 
erned by the strong interactions. In this paper we combine both methods follow- 
ing the pioneering work of Leinweber Utilising known results in QCD Sum 
Rules we apply them to Lattice Gauge Theory in order to check that the two are 
compatible. 

We begin by following the analysis of Refs. In this work, the Continuum 

Model (often used in QCDSR analyses p|) is tested using lattice Monte Carlo 
data. The Continuum Model replaces the discrete spectrum of excited states 
with a continuum of states of a certain spectral density. This spectral density is 
calculated using a Wilson OPE expansion of the relevant correlators in the small 
and intermediate time separation regime. It has been argued that this Continuum 
Model may be a good approximation to the actual situation of a large set of 
discrete excited states 

We extend the work of Refs. by including mesonic states. We also apply our 
analysis to several (3 values (corresponding to different lattice spacings), and to 
two actions with different lattice artefacts (namely the Wilson action and the 
"clover" Sheikholeslami-Wohlert (SW) improved action |^). Thus we are able to 
study the effects of lattice artefacts in the calculation. The lattice data we use 
in this study is from the APE Collaboration. As in the analysis of |]T|, we use 
zero-momentum correlation functions. At first sight the analysis suggests that 
the QCDSR Continuum Model is in reasonably good agreement with the lattice 
data, confirming the results of [|l|. However, this agreement is ambiguous since 
we also find that the naive (non-relativistic) quark model predicts essentially the 
same behaviour as the Continuum Model. Furthermore, because we have data at 
different lattice spacings, we are able to test the continuum scaling behaviour. 
We find that the parameter introduced to take account of lattice distortions in 
the Continuum Model does not scale towards its predicted value as the lattice 
spacing is reduced. Thus we conclude from our analysis that there are reservations 
regarding the validity of the QCDSR Continuum Model (but confirmation of the 
quark model). We also test the relativistic quark model against our data and find 
poorer agreement. 

Through this study of the QCDSR Continuum Model, we have defined a proce- 
dure which enables lattice correlation functions to be fitted at small (Euclidean) 
time values. This can be usefully employed in the case of correlation functions 
that have numerical problems in the isolation of the ground state properties at 
large time separations. By fitting the lattice correlation function data to a term 
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representing the ground state contribution, plus a term representing the contri- 
butions of the excited states, we show that accurate ground state parameters 
can be extracted even when the fitting window is very close to the source. We 
explicitly show that these QCDSR-inspired fitting functions reproduce the pa- 
rameters of the ground state better than the conventional, double-exponential 
fitting procedure. 

Recently, some other researchers have also used Continuum Model ideas to study 
lattice correlation functions They study hadronic correlation function, but 
use configuration space (i.e. point-to-point functions; see also f^) rather than the 
momentum space correlation functions analysed here. 

This paper is organised as follows: in Sect. 2 details of the lattice data used in 
our analysis are given; Sect. 3 describes the QCD Sum Rules approach; in Sect. 4 
results are presented for these fits; in Sect. 5 we present results using the quark 
model; and finally, in Sect. 6 there is a discussion and conclusion. An earlier 
version of this work appeared in 



2 Lattice Details 



The lattice data used in this study is from the APE collaboration and uses the 
quenched approximation. Specifically it is data which has been used primarily 



in the study of weak matrix elements such as fs and Bk Ww vM, and, more 



recently, for light hadronic spectroscopy . A list of the simulation parameters 
is presented in Table |l|. For full details of the simulations, see the corresponding 
reference. All errors in the results were obtained using the Jacknife method with 
10 configurations eliminated per cluster. 

The lattices in this study span a range of /5- values which correspond to different 
lattice spacings, a. In Table |I| and throughout this paper, a is set from the K 



and K* masses following the approach described in [jT3[. This method avoids any 
chiral extrapolation and uses data at the simulated values of the quark masses, 
i.e. at around the strange mass. It is therefore preferable than other methods, e.g. 
obtaining from the p-mass etc. Note that the lattice spacings are smaller than 
those used in Refs. [0 and and so clearly lattice artefacts should be smaller. 

Typically each simulation was performed with three different values of the hop- 
ping parameter, k. Only one of these (for each simulation) was chosen to be 
included in this study. This choice was made so that the vector meson masses 
in physical units from each simulation were equal to 1.1 GeV within errors (see 
Table |1|). The one exception to this is the Wilson case at /? = 6.1 where there 
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Table 1 

Simulation parameters used. The inverse lattice spacing a^^ was obtained from the K 
and K* masses [T^- 



Ref. P Action Number of Lattice a ^ a 

Configs Volume [GeV] [fm] 



Mv 
[GeV] 



llj 6.0 Clover 
lH 6.2 Clover 
|ll 6.4 Clover 



200 
200 
400 



18^ X 64 2.01(6) 0.098(3) 0.1425 1.10(3) 
18^ X 64 3.0(3) 0.065(6) 0.14144 1.14(11) 
24^ X 64 3.9(2) 0.050(2) 0.1403 1.12(4) 



IJ 6.0 Wilson 
|11| 6.1 Wilson 
111 6.4 Wilson 



200 18^ X 64 2.27(4) 0.087(2) 0.1530 1.11(2) 
230 18^ X 64 2.7(2) 0.073(4) 0.1510 1.27(8) 
400 24^ X 64 4.1(2) 0.048(2) 0.1492 1.10(3) 



was no K value available to match this requirement. 

The (zero-momentum) hadronic correlation function is defined for mesons as 
G^it) ^Y.{0\T{Jix,t)J{6,0m. (1) 



While the definition of the meson correlators is standard, the definition of the 
nucleon and delta correlation functions is specific to the APE simulations, and is 
given by 



1+74 



(0|T{J(f,t)J(0,0)}|0) 



(2) 



where the projector (1 + 74) onto positive parity states is introduced. Another 
common definition for the nucleon correlation function is the one used by the 
UKQCD Collaboration (see for example [0,0]). It is given by 



(0|T{J(f,t)J(0,0)}7MX^|0) 



(3) 



Of specific interest to this work is the definitions of the interpolating operators, J, 
defined below. The channels analysed in this work are the nucleon, pseudoscalar 
meson, (spatial component of the) vector meson and (temporal component of 
the) axial meson. The other channels are included below for completeness. 



4 



Nucleon: 



J(f, t) = , t)Cj5d\x, t)]u's{x, t)eabc (4) 

<5 



Delta: see 15 



3 

J(f, t) = E E[«"(^' t)C-fy{x, tMix, t)eabc (5) 
Scalar meson: 

J{x,t) =lf{x,t)d''{x,t) (6) 
Pseudoscalar meson: 

J(f, t) = Tt{x, t)75C?"(f , t) (7) 
Vector meson (spatial components): 

J{x,t) = Wv^{x,thd''{S,t) (8) 

Vector meson (temporal component): 

J(f, t) = Tt{x, t)-fid''{x, t) (9) 
Axial meson (spatial components): 

AS,t) = i^M"(f,t)757.rf"(^,t) (10) 

i 

Axial meson (temporal component): 

J(f,t) =M"(f,t)7574rf''(x,t). (11) 

In these equations, the indices a, b, c refer to colour, 5 is a spinorial index and 
C = 7472 is the charge conjugation matrix. 

All the baryonic correlation functions presented in this work have the sum over x 
in Eq. (§) replaced by a sum over every 3rd lattice site in each of the three spatial 
directions. This procedure, called "thinning", is due to the limitations imposed 
on the calculation by the size of the APE memory. It introduces extra states in 
the correlation functions which are not present when the usual full sum over x 
it taken. These unwanted states have large energy corresponding to large spatial 
momentums and do not affect the conventional extraction of the ground state 
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properties [|T6|,|T0|,|T3[ . This is because at the (relatively large) times where there is 
ground state dominance, the higher energy states have decayed away. However, 
since the procedure presented in this paper uses small time values, thinning does 



have an effect which will be referred to in Sect. 4.3 



3 QCD (Sum Rules) Continuum Model 



In this section we describe the QCDSR Continuum Model. The basic object in this 
approach is the quark propagator in Euclidean space whose first few polynomial 
terms in the Wilson OPE expansion read (in the coordinate gauge, x'^A^^ = 
0), 



^ 27r2 x4 47r2 8tt^ 

-^{■.qq:)6'''^' + ■■■. (12) 

From Eqs. (0,0 & ^, the timesliced (three- momentum-projected) two-point cor- 
relation function then has the OPE expansion 

"0 1 

Gr''W= E -■aO„(m„(:gg:)), (13) 

n=— oo 



where C„ is a numerical coefficient, 0„ is some function, and no a positive integer. 
For the nucleon and delta, this leading term turns out to be no = 6 for the 
APE definition Eqs. (^. However, different definitions of the correlators lead to 
different QCDSR Continuum Model formulae, and if we were to use the UKQCD 
definition, Eq. (^, we would have no = 5 for both the nucleon and delta. For the 
mesons the exponents of the leading terms are given by no < 3. 

It is now useful to express the timesliced correlation function, 6*2(^)5 in the spec- 
tral representation 

00 

G2{t) = Jp{s)e-'^Us. (14) 
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The OPE expansion of the spectral density function, p'-^^^{s), is calculated by 
means of a simple inverse Laplace transform of Eq. (0) , 



pOPE^^^ = E r^TTj ■ CnOnim,, (: qq :)), (15) 

n=i 



where Re{t) > 0. In the sum over n, we have kept only the positive values of n 
since we are interested in the leading terms as t — > 0. 

The QCD Continuum Model is introduced at this stage by setting a threshold Sq 
in the energy scale s in Eq. (|Hp, so that the excited states' contribution to G2(t) 
is given only by the energies above sq. Performing the integral over s we obtain, 



GT*(t) = I p^^^(s) e-^* ds 



so 

no n-1 1 k 

= E E Zyi^ • CnO.(m„ (: qq :)). (16) 

n=l k=0 

Thus the continuum contribution to the correlation function has been derived us- 
ing QCDSR style considerations. The hope is that it correctly models the discrete 
states in the spectrum above the ground state. 

However, the full correlation function contains the ground state as well as the 
above continuum contribution Gf"'*. In the quenched approximation, the ground 
state contributes a delta function to the full spectral density ||l|J 



Pis) = ^S{s - M) + e{s - .o)p^^^(^) (17) 



where p^^^{s) is given by Eq. (p!5|). The full correlation function can then be 
written as (from Eqs. (|l^ & |16D) 



Z 



G2(t) = ^e-*^*+e/p^^^(.)e-- d. 



so 



y no n-1 i k 

e-^*+ee--* E ■ C^nO„(m„ (: qq :)) 



2M k\ t 



There are four parameters in this Ansatz for G2{t): 
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Z is the normalization of the ground state. Note that in |T]J^ the parameter Ai 
was used. Ai and Z are related by Af = Zl{2M)] 
M is the mass of the ground state particle; 

^ is a new parameter introduced to normalise the contribution of the excited 
states n. In the continuum limit, ^ should be equal to one. (We will discuss 



this in more detail in Sect. [4.3|) ; 
- So is the continuum threshold which parametrises the onset of the excited 
states. 

We have developed some FORM codes able to calculate the coefficients C„ in the 
expansions in Eqs. (^) and (|18|) for the operators in Eqs. - |1T]). The relevant 
OPE expansions (i.e. the G^^^{t) in Eq. (|13D) turn out to be: 



Nucleon (APE data): 

2 ^ ' 2567r4 16 647r4 16%^ 
7 (: qq :) 



327r2 t3 
Delta (APE data): 



+ ■■■ (19) 



15 (: qq :) 



47r2 t3 

Scalar meson: 



+ ■ ■ • (20) 



Pseudoscalar meson: 

Vector meson (Spatial components): 



Axial meson (Temporal components): 

—3 m 
2^1 



Gr^W = 7r4-^ + --- (24) 
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Axial meson (Spatial components): 



qOPE 



it) 



1 1 



3 mi 



2ttH^ 27r2 t 



+ 



(25) 



Note that for the vector meson, the leading term for the temporal component is 
0{m^). Since this corresponds to the first neglected term in the OPE expansion 
for the quark propagator in eq. ([l2|) , it is beyond the scope of this work and we 
therefore do not consider this channel again. Note also that there is no term 
0{mg) for the spatial component of the vector correlation function. 

In the nucleon and delta case we have included for completeness the contribution 
of the vacuum condensate. However, in practice, it is not necessary to consider 
neither this term nor the m^/t^ term, since they are numerically irrelevant. We 
note that Leinweber came to the same conclusion . It is important to note that 
for all mesons listed there is no term appearing which is one order in t smaller 
than the leading term. i.e. the first nonzero correction, with respect to the leading 
term, is of order and not t as we would naively expect. Similarly, the order 
t^ term (with respect to the leading term) is zero. The first appearance of the 
vacuum condensate (: gg :) in the meson case is of order relative to the leading 
term, rendering its contribution totally insignificant. 

We give here for the sake of completeness the correlation functions corresponding 
to Eq. (p!^), i.e. the G2{t) actually used to fit the lattice data. They are: 



Nucleon (APE data): 

Z 



G2it) 



Mt 



2M 

+^e-'o* 




1 s 



1 So Isl 

- + — H H H + 

t^ 2 1^ Q t^ 2At^ 120 t 







t^ t^ 2t^ Gt'' 24:t 



1 So 

h — + 

t^ t^ 2t^ 




(26) 



Delta (APE data): 

G2{t)- ^ 




1 



t^ 2 1^ Qt^ 2At^ 120 t 



1 Sn 1 sfi 1 Sn 

h — H H + 

t^ t^ 2t^ Qt^ 2A t 
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Scalar mesons: 



G2it) 



Z 



-Mt 



2M 

Pseudoscalar mesons 

Z 



-sot 



G2{t) 



2M 



■sot 



AnAt^ 2 t , 



AnAf^ 2 t 



Vector mesons (Spatial components): 



-Mt 



2M 



-so* 



1 



1 So ISg' 

2^lt3 + + 27 



Axial mesons (Temporal components): 



-Mt 



2M 



-sot 



3 

2^T 



Axial mesons (Spatial components): 



G2{t) 



-Mt 



2M 



-sot 



27r2 \ t3 



So 



IsJ 
2 t 
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4^T 



3 

4^T 



3 

2^T 



(27) 



(28) 



(29) 



(30) 



(31) 



We have carried out several four-parameter fits to the nucleon, pseudoscalar, vec- 
tor (spatial components) and axial-vector (temporal components) using Eqs. ( p6| , 
[28| , |29| & |3^) to check the ideas presented in this paper. The results of these fits 
are described in detail in the following Section. 

One could also introduce a further parameter A as upper limit of integration over 
s in Eq. (|T8|) to account for the lattice cut-off. Ref. ^ reports that this does not 
make any difference as long as t is not taken less than two. For this reason we 
restrict our fits to i > 2. 



4 Fits to data 



Correlation functions from the simulations outlined in table ^ were fitted to the 
following three different functional forms: 

Conventional Single State Fit ("1-exp") 

F{t) = ^e-"' (32) 
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This is the usual functional form used in the study of lattice correlation func- 
tions. 

QCD Continuum Model Fit ("Cont") 



The fitting forms in Eqs. ( P6| - pl| ). (i.e. those obtained from the O.P.E. plus the 
Continuum Model assumption in Sect. ^) 
Conventional Two State Fit ("2-exp") 

F(t) = —e-''' + ^e-^'* (33) 

This functional form was chosen since it is traditionally used as the generali- 
sation of Eq. ( P^ ) when attempts are made to include the effects of the higher 
mass state(s). It is an alternative to the Continuum Model fits (i.e. it has 
the same number of parameters) and the results of both fits will be directly 
compared. 

Note that for the mesons, we have symmetrised the ground state exponential in 
time for all three cases above by taking into account the backward moving state 
in the fits. 



The conventional single state fits are only valid for asymptotic states. This limits 
their applicability to a time window where the "effective mass" is constant in 
time.i We use the results of these fits (in the asymptotic region) as standard 
values of the ground state fitting parameters Z and M with which to compare 
the results of the other two fitting methods. A list of the parameters Z and M and 
the obtained from these single state fits, is displayed in Tables @,^|| & |^ in the 
rows marked "1-exp". The four channels: nucleon and the pseudoscalar, vector 
(spatial components) and axial (temporal component) mesons were analysed. 

The time windows used in each of the fitting methods are listed in Table |^ and are 
chosen so that the effective masses are stable. The starting points in the windows 
are a function of the (3 value, as expected, since the lattice spacing varies with /3. 
In the case of the QCD Continuum Model fits, we use time windows which begin 
very close to the origin: t = 2 - 28 for /5 = 6.0; t = 3 - 28 for /5 = 6.1 & 6.2; 
t = 4 — 28 for j3 = 6.4. (See Table ^.) The starting times for these windows are 
chosen with the criterion that t >2 (see Sect. ^) and that they are approximately 
equal in physical units for all simulations. The results of these fits are displayed 
in Tables & ^ in the rows marked "Cont". 

So that a direct comparison can be performed, the same time windows were used 
to perform conventional two state fits using Eq. (P3|). The results of the fitting 
parameters from these fits are also displayed in Tables || & ^ in the rows 

^ The effective mass is defined as log{G2{t)/G2{t + 1)). 
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labelled "2-exp". 



We now discuss the results of these fitting methods in detail. 
4-1 Ground State Parameters 

Concentrating on the comparison of the ground state parameters, Z and M, from 
the "Cont" fit with those from the "1-exp" fit, we see good agreement for the 
mesonic states. This is particularly encouraging because very small time values 
were included in the fitting window. However, for the baryons, poorer agree- 
ment is observed. In fact the agreement between the Z and M values from the 
Continuum Model and "1-exp" fit becomes worse as (3 increases (i.e. as the con- 
tinuum limit is approached). Overall, these fits suggest that the QCD Continuum 
Model is a reliable method of extracting ground state information in the mesonic 
case, but that there is some question about its validity for the baryons. Note that 
Leinweber, who studied baryons only, found good agreement flj. A possible ex- 
planation for this is that he used data at a larger lattice spacing and with poorer 
statistics than those studied here. Note also that the baryonic data is "thinned" 
which will obviously affect the normalisation of the higher mass states. This could 
conceivably also affect the ground state fitting parameters in the "Cont" fits. 

The results for Z and M for the "2-exp" fits require careful interpretation. Com- 
paring the "2-exp" and "1-exp" fits we see that the values for both Z and M 
for the "2-exp" fits are larger than those in the "1-exp" fit in every case. This 
presumably means that the introduction of a "raw" second exponential biases 
the parameters of the asymptotic state. Also, it can be seen that in every case 
except one (the Wilson (3 = 6.0, Pseudoscalar meson) the "Cont" values for Z 
and M reproduce the "standard values" (i.e. the "1-exp" values) better than the 
"2-exp" do. 

We can understand this bias in the case of the "2-exp" fit as follows. The true set 
of states that contributes to the spectrum includes the ground state and many 
excited states. In the case of the "2-exp" fit, all the excited states are modelled as 
a single state. Clearly the values of M' and Z' from this fit will be some average 
over all the excited states. This means that the M' value will be larger than the 
first excited state's mass. However, of the excited states, it is this first excited 
state that is the most significant correction to the ground state at intermediate 
values of t (i.e. at t values large enough so that the excited states are no longer 
dominant, and small enough so that the ground state also has not yet dominated 

^ Note also that it is well known that the statistical errors in the fitted Z and M 
values typically underestimate the real errors due to correlations in the data |p^]. 
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the correlation function). This means that the correlation function for the "2-exp" 
fit will be too small in these intermediate values of t. To adjust for this mismatch, 
the value of the ground state's mass, M (and Z since they are correlated) will be 
shifted slightly higher in the "2-exp" fit. 

Obviously, in order to avoid this bias, correct account must be taken of all the 
excited states in the spectrum. The Continuum Model is one way of modelling 
these states, and thus should overcome at least part of the above shortcomings 
of the "2-exp" fit. The agreement of the ground state parameters, Z and M, 
between the Continuum Model and "1-exp" fit confirms this hypothesis. 



4-2 Quality of the Fits 

Turning to the values, the first thing to note is that the x^/ d.o.f. for the "Cont" 
and "2-exp" fits can be as much as 10 or 100. At first sight this is a worrying 
feature, especially since the is calculated ignoring correlations between times- 
lices. (The presence of correlations canonically lowers the uncorrelated x^/d-o.f. 
to below unity.) However, the correlation function data at small times has rela- 
tively tiny errors compared with large times. (Typically the correlation function 
has statistical errors of one part in a thousand at the beginning of the fitting 
window, and around one part in ten at the end of the window.) Therefore any 
small discrepancy between the fitting functions and the data at small times leads 
to a very large contribution in the x^- It is unreasonable to expect that the small 
time behaviour of the correlation functions could possibly be reproduced at the 
level of these tiny statistics by any model containing only four parameters, no 
matter how well physically motivated. It is also worth bearing in mind that the 
meson correlation functions fall by from 3 to as much as 7 orders of magnitude 
in the region fitted in the "Cont" and "2-exp" cases. In the case of the baryons 
the correlation functions can fall by as much as 13 orders of magnitude. Fig. |I] 
shows the (natural) logarithm of the correlation function for the vector meson in 
the Wilson, j3 = 6.4 case (which is a representative mesonic example). The plot 
shows the fits from all three methods above. This figure illustrates the points 
made above; the errors in the data for small times are clearly tiny. 

We have argued that the absolute values of the x^/d.o.f. values might reasonably 
be expected to be large for the "Cont" and "2-exp" fits. However their relative 
values should correspond to the relative quality of fits of the two methods. From 
Tables & |^ we see that the for the "Cont" fit are smaller than (or, at 
worse, similar to) the "2-exp" fit in every case. In some cases (particularly for the 
Pseudoscalar meson) the value for the "Cont" case is as much as an order of 
magnitude smaller than the "2-exp" case. This is again evidence that the "Cont" 
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fitting functions are an improvement over tlie traditional "2-exp" functions. 



As examples of the quality of the fits, Figs. ^ |^ and ^ show the effective mass 
plots for the Clover, /? = 6.2, Pseudoscalar case, the Wilson, /? = 6.4, Vector 
meson case, and the Wilson, /? = 6.0, nucleon case respectively. The central 
values of the effective mass from the three methods: "1-exp"; "Cont"; and "2- 
exp" are shown. (Note that the effective mass has been corrected for the effect 
of the backward moving state and should therefore be constant at large times 
approaching the middle of the lattice.) In these figures the "2-exp" fit can be seen 
to give a reasonable reproduction of the Monte Carlo data, however, the "Cont" 
method performs a little better. In particular the effect discussed in Sect. at 
intermediate time values can be seen in the "2-exp" case. These three cases are a 
representative sample of all the cases studied, and support the conclusions above 
regarding the relative merits of the "2-exp" and "Cont" approaches. 



Further evidence that the "Cont" method is an improvement over the "2-exp" 
case can be seen in Figs. ^ and ^. The channels shown are for the pseudoscalar 
and nucleon respectively for the Clover, (3 = 6.0 simulation. In these figures the 
mass, obtained from both the "Cont" and "2-exp" fits using t = 2 — > tMAx, is 
plotted against Imax- (The "1-exp" fit using t = 12 — 28 is shown as the solid 
horizontal line.) As can be seen from the fits, the "Cont" method converges much 
faster to the true mass value compared with the "2-exp" case. Other channels 
show a similar behaviour. 



4-3 Parameters of the Excited State(s) 



Clearly it is not sufficient for the "Cont" fits to reproduce the ground state 
parameters Z and M, and to have a sensible x^, they must also give values for 
the continuum parameters, Sq and ^, which are acceptable within the assumptions 
of the QCDSR Continuum Model. 

In Table |^ the values for sq (in GeV) are listed for all the cases studied. These 
values must satisfy two criteria: (i) since Sq corresponds to a physical threshold, 
they should be constant (in GeV) for each channel as P is varied; and (ii) they 
should be large enough to be in a region where perturbation theory is valid. The 
second criterion is required so that the perturbative expression for the quark 
propagator, Eq. (0), is valid. From Table |^, it is clear that the Sq values are 
roughly constant in (3 for the mesons, but not for the nucleon (particularly the 
Clover case). Also, for the pseudoscalar meson, the Sq values are arguably too 
small ( ^ 2 GeV) for perturbation theory to be considered reliable. Thus only 
the vector and axial mesons convincingly pass the criteria above. 
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We now turn to the values for ^. Note that this parameter was introduced to 
parametrise the distortions introduced in the lattice formalism. It should therefore 
be unity in the continuum limit Table || displays the various values of ^ 
from the "Cont" fits. Clearly, in the case of the axial meson, the ^ values are 
not "~ 1 + (9(a)". The large size of these ^ is because the leading term in the 
Continuum part of G2{t) for this channel is 0{w?q) (see Eqs. (0) and (|2^)). This 
implies that ^ values several orders of magnitude larger than the other channels 
(which have 0{1) terms in G2{t)) are required. 

In conclusion, the results of the continuum parameters, Sq and ^, cast doubt over 
the applicability of the "Cont" model, since the values obtained from the fits are 
not consistent with the assumptions made in the QCDSR Continuum Model. 

Note, it is not possible to compare directly the values for ^ for Wilson, nucleon 
data with those in Ref. This is because the baryonic channels studied here 
have "thinned" correlation functions (see the discussion in Sect. This affects 
the short-time region only, due to the presence of extra high energy states, and 
thus conceivably alters the values of ^. 



4-4 Discussion on the Fits 



Summarising the subsections above, we see that the QCDSR Continuum Model 
fitting functions reproduce the correct ground state parameters, Z and M, of the 
Monte Carlo data in the meson case. The Continuum Model is not as successful 
however in the baryonic sector. The "2-exp" case does not appear to reproduce 
the ground state parameters as effectively as the "Cont" case for both baryons 
and mesons. We presented an explanation for this finding in Sect. [4.1| . However, 
the values found for the continuum parameters Sq and ^ cast doubt over the 
applicability of the "Cont" Ansatz, since the values obtained are not consistent 
with the assumptions used in its derivation. The conclusion therefore is to reject 
the QCDSR Continuum Model as a way of modelling these lattice correlation 
functions. 



However it is clear that, apart from taking on non-physical values of its fitting 
parameters, the "Cont" Ansatz does fairly adequately reproduce the lattice cor- 
relation function data. Therefore we require similar fitting Ansatze as those in 
the "Cont" case, but without the corresponding restrictions on the parameters 
So and ^. 

As we shall see in the next section, the naive non-relativistic quark model predicts 
essentially the same functional form for the hadronic correlators as in the QCDSR 
Continuum Model (for both the mesonic and baryonic cases). Therefore it is a 
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15 

t 



Fig. 1. Plot of the (natural) logarithm of the correlation function for the (3 = 6.4, 
Wilson, Vector meson data. The curve is the result of the "Cont" fit and the straight 
line the result of the "1-exp" fit. The "2-exp" fit is indistinguishable from the "Cont" 
fit on this plot. 



candidate model to use to fit the data. 

Before commencing tlie discussion on tlie quark model, we comment on the results 
of the Wilson and clover actions. It is natural to ask if we can observe any 
difference in the fitted parameters for these two actions which might be a signal 
for improvement. Studying the values Tables ^ & |^ it is difficult to 

see discernible differences between the two actionsTa There are, however, clear 
differences between the Continuum parameters in the two actions, particularly 
for ^ (see Table This is not surprising since differences between the two actions 
are known to occur at small times values (see 0). However, the interpretation of 
these effects will require more work. 



It is possible that there is some difference in tlie vector channel - see Table |^. 



16 



Nucleon 



Clover Action 



f3 




6.0 


6.2 


6.4 


Z 


1-exp 


[0.14 ± 0.02] X 10-5 


[0.14 ±0.05] X 10-6 


[0.15 ±0.04] X 10- 




Cont 


[0.24 ± 0.02] X 10-5 


[0.30 ± 0.03] X 10-*^ 


[0.48 ± 0.02] X 10" 




2-exp 


[0.34 ± 0.04] X 10-5 


[0.33 ± 0.03] X 10-*^ 


[0.50 ± 0.02] X 10- 


Ma 


1-exp 


0.828 ± 0.009 


0.59 ±0.02 


0.438 ± 0.009 




Cont 


0.865 ± 0.008 


0.636 ± 0.009 


0.491 ± 0.003 




2-exp 


0.89 ±0.01 


0.643 ± 0.009 


0.492 ± 0.003 




1-exp 


(l.±l.)/15 


(0.7±0.8)/9 


(0.02±0.05)/3 




Cont 


(1000. ± 70.)/23 


(310. ±60.)/22 


(1000. ±60.)/21 




2-exp 


(4000. ± 100.)/23 


(550. ±80.)/22 


(1000. ± 70.)/21 



Wilson Action 



/3 




6.0 


6.1 


6.4 


Z 


1-exp 


[0.16 ± 0.02] X 10-5 


[0.61 ±0.07] X 10-6 


[0.23 ± 0.04] X lO-'^ 




Cont 


[0.24 ± 0.02] X 10-5 


[0.127 ±0.008] X 10-5 


[0.92 ±0.05] X lO-'^ 




2-exp 


[0.32 ± 0.03] X 10-5 


[0.144 ±0.009] X 10-5 


[0.105 ±0.006] X 10-' 


Ma 


1-exp 


0.797 ± 0.008 


0.737 ± 0.006 


0.428 ± 0.006 




Cont 


0.823 ± 0.007 


0.776 ± 0.005 


0.495 ± 0.004 




2-exp 


0.845 ± 0.008 


0.784 ± 0.005 


0.503 ± 0.005 




1-exp 


(0.3 ±0.5)/15 


(0.1 ±0.3)/ll 


(0.02 ±0.03)/3 




Cont 


(220. ±40.)/23 


(660. ±60.)/22 


(2000. ±90.)/21 




2-exp 


(1000. ±80.)/23 


(1000. ±80.)/22 


(2000. ± 100.)/21 



Table 2 

Values for the fitting parameters for the nucleon. The 1-exp, Cont and 2-exp refer to 
fits using Eqs. p^, and respectively. 
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Pseudoscalar Meson 



Clover Action 



/3 




6.0 


6.2 


6.4 


Z 


1-exp 


0.0328 ± 0.0009 


[0.80 ± 0.05] X 10-2 


[0.27 ±0.01] X 10-2 




Cont 


0.0329 ± 0.0007 


[0.86 ± 0.03] X 10-2 


[0.278 ± 0.009] X 10-2 




2-exp 


0.0349 ± 0.0007 


[0.92 ± 0.03] X 10-2 


[0.312 ± 0.008] X 10-2 


Ma 


1-exp 


0.438 ± 0.001 


0.294 ± 0.003 


0.220 ±0.001 




Cont 


0.438 ± 0.001 


0.296 ± 0.002 


0.220 ±0.001 




2-exp 


0.441 ± 0.001 


0.299 ± 0.002 


0.224 ± 0.001 




1-exp 


(0.2±0.3)/15 


(0.01 ±0.05)/9 


(0.001 ±0.003)/3 




Cont 


(3. ± l.)/23 


(3. ±2.)/22 


(1.1 ±0.5)/21 




2-exp 


(48. ±5.)/23 


(18. ±4.)/22 


(21. ±4.)/21 



Wilson Action 







6.0 


6.1 


6.4 


z 


1-exp 


0.0244 ± 0.0007 


0.0149 ± 0.0005 


[0.210 ± 0.009] X 10-2 




Cont 


0.0226 ± 0.0007 


0.0144 ± 0.0004 


[0.200 ± 0.009] X 10-2 




2-exp 


0.0253 ± 0.0006 


0.0155 ± 0.0004 


[0.236 ± 0.005] X 10-2 


Ma 


1-exp 


0.423 ±0.001 


0.404 ± 0.001 


0.205 ± 0.001 




Cont 


0.420 ± 0.001 


0.403 ± 0.001 


0.204 ± 0.002 




2-exp 


0.425 ± 0.001 


0.406 ± 0.001 


0.210 ±0.001 




1-exp 


(0.3±0.4)/15 


(0.1 ±0.2)/ll 


(0.002 ±0.004)/3 




Cont 


(20. ±4.)/23 


(1. ± l.)/22 


(0.4±0.6)/21 




2-exp 


(23. ±4.)/23 


(18. ±5.)/22 


(27. ±5.)/21 



Table 3 

Values for the fitting parameters for the pseudoscalar meson. The 1-exp, Cont and 2-exp 



refer to fits using Eqs. (3i), ^q) and i\33l) respectively. 
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Vector Meson (Spatial Components) 



Clover Action 







6.0 


6.2 


6.4 


z 


1-exp 


[0.93 ± 0.06] X 10-2 


[0.18 ±0.03] X 10-2 


[0.51 ±0.04] X 10 




Cont 


0.0106 ± 0.0003 


[0.27 ±0.01] X 10-2 


[0.77 ± 0.03] X 10 




2-exp 


0.0121 ± 0.0004 


[0.300 ± 0.010] X 10-2 


[0.94 ± 0.03] X 10 


Ma 


1-exp 


0.548 ± 0.004 


0.378 ± 0.006 


0.285 ± 0.003 




Cont 


0.556 ± 0.003 


0.397 ± 0.003 


0.300 ± 0.002 




2-exp 


0.565 ± 0.003 


0.403 ± 0.003 


0.309 ± 0.002 




1-exp 


(l.±l.)/15 


(0.4±0.6)/9 


(0.02±0.02)/3 




Cont 


(130. ±30.)/23 


(120. ±20.)/22 


(91. ±2.)/21 




2-exp 


(610. ±60.)/23 


(270. ±40.)/22 


(290. ±30.)/21 



Wilson Action 



6.0 



6.1 



6.4 



Z 1-exp 0.0102 ± 0.0004 

Cont [0.95 ± 0.03] x 10-2 

2-exp 0.0120 ± 0.0003 

Ma 1-exp 0.508 ± 0.003 

Cont 0.505 ± 0.002 

2-exp 0.517 ± 0.002 

x2 1-exp (0.3±0.6)/15 

Cont (20. ±5.)/23 

2-exp (200. ±20.)/23 



[0.57 ±0.03] X 10-2 
[0.59 ±0.02] X 10-2 
[0.68 ± 0.02] X 10-2 
0.466 ± 0.002 
0.468 ± 0.002 
0.474 ± 0.002 
(0.4±0.6)/ll 
(8. ± 5.)/22 
(110. ±20.)/22 



[0.62 ± 0.04] X 10-3 
[0.74 ± 0.03] X 10-3 
[0.99 ± 0.03] X 10-3 
0.267 ± 0.003 
0.274 ± 0.002 
0.286 ± 0.002 
(0.008 ±0.009)/3 
(11. ±4.)/21 
(180. ±20.)/21 



Table 4 

Values for the fitting parameters for the (spatial components of the) vector meson. The 



1-exp, Cont and 2-exp refer to fits using Eqs. plj, (2i) and respectively. 



19 



Axial- Vector Meson (Temporal Components) 



Clover Action 



/3 




6.0 


6.2 


6.4 




z 


1-exp 


[0.171 ± 0.008] X 10-2 


[0.37 ±0.03] X 10-3 


[0.111 ±0.009] X 10~ 






Cont 


[0.179 ±0.006] X 10-2 


[0.38 ± 0.02] X 10-3 


[0.118 ±0.004] X 10- 


-3 




2-exp 


[0.180 ±0.006] X 10-2 


[0.38 ± 0.02] X 10-3 


[0.118 ±0.004] X 10" 


-3 


Ma 


1-exp 


0.437 ±0.002 


0.293 ± 0.004 


0.217 ±0.003 






Cont 


0.440 ± 0.002 


0.295 ± 0.004 


0.219 ±0.002 






2-exp 


0.440 ± 0.002 


0.295 ± 0.004 


0.219 ±0.002 






1-exp 


(0.3±0.2)/15 


(0.06±0.06)/9 


(0.02 ±0.02)/3 






Cont 


(10. ±3.)/23 


(1.2±0.8)/22 


(2. ± l.)/21 






2-exp 


(15. ±4.)/23 


(1.5±0.9)/22 


(2. ± l.)/21 





Wilson Action 



f3 




6.0 


6.1 


6.4 


Z 


1-exp 


[0.23 ± 0.01] X 10-2 


[0.159 ±0.006] X 10-2 


[0.126 ±0.009] X 10-3 




Cont 


[0.234 ± 0.007] X 10-2 


[0.170 ±0.004] X 10-2 


[0.140 ±0.005] X 10-3 




2-exp 


[0.236 ± 0.006] X 10-2 


[0.171 ±0.004] X 10-2 


[0.140 ±0.005] X 10-3 


Ma 


1-exp 


0.422 ±0.002 


0.403 ± 0.002 


0.201 ± 0.002 




Cont 


0.423 ±0.002 


0.406 ± 0.002 


0.205 ± 0.001 




2-exp 


0.423 ± 0.002 


0.406 ± 0.002 


0.205 ± 0.001 




1-exp 


(0.5±0.5)/15 


(0.8±0.5)/ll 


(0.01 ±0.01)/3 




Cont 


(2. ± 2.)/23 


(5. ±3.)/22 


(0.8±0.9)/21 




2-exp 


(4. ± 2.)/23 


(6. ±3.)/22 


(0.9± 1.0)/21 



Table 5 

Values for the fitting parameters for the (temporal components 



1-exp, Cont and 2-exp refer to fits using Eqs. ^M), (3C) and 



of the) axial meson. The 
respectively. 
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Fig. 2. Effective mass plot for the (3 = 6.2, Clover, Pseudoscalar data. The values of 
M from the fits to the "1-exp" is shown by a horizontal solid line. The "Cont" is the 
dashed curve, and the "2-exp" is the dotted curve. 




Fig. 3. Effective mass plot for the (3 = 6.4, Wilson, Vector meson data. The values of 
M from the fits to the "1-exp" is shown by a horizontal solid line. The "Cont" is the 
dashed curve, and the "2-exp" is the dotted curve. 
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Fig. 4. Effective mass plot for the f5 = 6.0, Wilson, Nucleon data. The values of M from 
the fits to the "1-exp" is shown by a horizontal solid line. The "Cont" is the dashed 
curve, and the "2-exp" is the dotted curve. 



0.52 




0.46 



Fig. 5. Results for the mass from the fit to the (3 = 6.0, Clover, pseudoscalar case for 
times i = 2 — > tMAX- The "Cont" fit results are shown as diamonds, and the "2-exp" 
case as squares. The "1-exp" fit for t = 12 — > 28 is depicted by a horizontal solid line. 
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Clover Action 



P 6.0 

1- exp 12-28 
Cont 2-28 

2- exp 2-28 



6.2 6.4 

18-28 24-28 

3-28 4-28 

3-28 4-28 



Wilson Action 



p 


6.0 


6.1 


6.4 


1-exp 


12-28 


16-28 


24-28 


Cont 


2-28 


3-28 


4-28 


2-exp 


2-28 


3-28 


4-28 



Table 6 

The time windows used in each of the fitting methods. 
5 Quark Models 



5.1 Non-relativistic Quark Model 



The quark model has had considerable success in reproducing many of the features 



of the hadronic spectrum [18|. It is thus natural to check its predictions for the 
hadronic correlators with the lattice Monte Carlo data. We present firstly the 
derivation of the hadronic correlators for the non-relativistic quark model, and 
then the relativistic case. 



The aim of this exercise is to determine the density of states function p{s) which 
defines the number of states per unit energy range. Taking the quarks as on- 
shell, there are 3 degrees of freedom available for each quark in the hadron: one 
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so [GeV] 



Clover Action 

P 6.0 6.2 6.4 

Nucleon 5.8 ± 0.2 8.5 ± 0.8 10.1 ±0.5 

Pseudoscalar Meson 1.87 ± 0.06 2.0 ± 0.2 1.6 ±0.1 

Vector Meson 2.57 ± 0.09 2.9 ± 0.2 2.4 ±0.1 
(Spatial Components) 

Axial- Vector Meson 4.1 ± 0.2 5.0 ± 0.4 5.3 ± 0.5 
(Temporal Components) 



Wilson Action 



(3 6.0 6.1 6.4 

Nucleon 5.0 ±0.1 5.9 ± 0.4 6.9 ± 0.3 

Pseudoscalar Meson 1.34 ±0.05 1.7 ±0.1 1.4 ±0.1 

Vector Meson 1.77 ±0.04 2.0 ±0.1 1.97 ±0.08 
(Spatial Components) 

Axial- Vector Meson 3.5 ± 0.1 3.5 ± 0.3 4.2 ± 0.5 
(Temporal Components) 

Table 7 

Values for the fitting parameter sq in GeV from the "Cont" fit. 
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Clover Action 

P 6.0 6.2 6.4 

Nucleon 72.4 ±0.5 49.7 ±0.6 22.6 ± 0.4 

Pseudoscalar Meson 9.49 ± 0.07 6.8 ±0.1 4.87 ± 0.07 

Vector Meson 9.41 ± 0.07 5.70 ± 0.10 2.85 ± 0.04 
(Spatial Components) 

Axial- Vector Meson [1.54 ± 0.06] x 10^ [1.3 ± 0.2] x 10^ [6. ± 2.] x 10^ 
(Temporal Components) 



Wilson Action 

/3 6.0 6.1 6.4 

Nucleon 6.73 ± 0.03 5.58 ± 0.06 1.82 ±0.04 

Pseudoscalar Meson 3.82 ± 0.01 4.03 ± 0.06 3.74 ± 0.04 

Vector Meson 3.40 ± 0.02 3.28 ± 0.04 2.66 ± 0.02 

(Spatial Components) 

Axial- Vector Meson 126. ± 8. 56. ± 2. 94. ± 8. 

(Temporal Components) 

Table 8 

Values for the fitting parameter ^ from the "Cant" fit. 
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1.6 



15 20 
t Max 



Fig. 6. As in Figure ^ but for the Nucleon case. 

per spatial momentum component. Therefore in the case of mesons there are 
naively 6 degrees of freedom, with 9 degrees of freedom for baryons. However, 
there are 4 constraint equations. Three correspond to demanding that the total 
3-momentum is fixed to the hadron's momentum (which is zero in the correlators 
we are studying), and one which corresponds to demanding that the total kinetic 
energy is fixed to s. This leaves 2 degrees of freedom for the mesons and 5 degrees 
of freedom for the baryons. This simple, non-relativistic analysis leads to the 
following form for the spectral density. 



For mesons: Pnrqm{s) ~ (34) 
For baryons: pnrqm{s) ~ (35) 

The above continuum-like behaviour for p{s) is clearly most appropriate for large 
s where the quarks are approximately free. Therefore the actual density of states 
can be approximated by 

P(^) = 7777 '^(S - ^) + ~ So)pnrqm{s), (36) 



where the delta-function represents the ground state, and the non-relativistic 
quark model result is used for the continuum (beginning at the threshold energy 
So). Using this definition of p{s) in the spectral equation (Eq. ([l^)) "we obtain 
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for mesons: 

^ A,-.. . A- . I 4) e-'. (37) 



and for baryons: 

CM - ^e-" + (i + I + A , A ^ _£L ^ A) ,38) 



where and K' are some numerical constants. Note that Eqs. (P7|fcp^) are iden- 
tical to Eqs. p6| - |3TD in Sect. ^ with = 0. Since numerically the terms 0{mg) 
and (9(m^) are insignificant, we have shown that the non-relativistic quark model 
predicts essentially the same functional form as the QCD Sum Rules Continuum 
Model for these channels. 



Thus we have achieved the aim raised in Sect. [4.4| : we have found an Ansatz 
which reproduces the data better than the "2-exp" fits, but which doesn't suffer 
from the unphysical parameter values of the "Cont" fit. 



5.2 Relativistic Quark Model 



We now turn our attention to the derivation of p{s) in the relativistic quark 
model. In this case we can represent the number of states as: 

Number of states ~ / ^ / ^^'(Pi +P2)5(rf^)^ + (pf)^ - s% (39) 



in the meson case, with an obvious generalisation for baryons. The main differ- 
ence between this case and the non-relativistic case is the energy factors in the 
denominators of the normalisation. Thus we obtain two powers of s less in the 
meson case, and 3 powers of s less in the baryonic case. i.e. 



For mesons: Prqm{s) ~ s (40) 
For baryons: Prqm{s) ~ (41) 

Following the above analysis, the 2-point hadronic correlators can easily be de- 
rived; 
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for mesons: 

G,{t) = Ae-A« + K Q e--* (42) 
and for baryons: 

G.it) ^ ^e- + /r (i + I + I) e-'K (43) 

Note that the relativistic quark model prediction for G2{t) for baryons is identical 
to the non-relativistic mesonic 6*2 (^)- 

5.3 Quark Model Fits 

It is natural to wonder if there is any way of using the lattice correlation function 
data to pin down the best form of the density of states p{s). With this in mind 
we assume that the density of states has the following form: 

^^'^ = ~ + ~ ^ ^^^^ 

where n is to be determined from the fit. The two-point correlation function for 
this Ansatz is defined, as usual, from eg. (p!^ . Eq. (^) has both the non-relativistic 
and relativistic forms as special cases. Fig. |^ shows the value plotted against 
n for fits to the 6 sets of pseudoscalar data in table |I| using eq.(|^). As can be 
seen from the plot, there is a distinct minimum at = 2 for five of the six data 
sets. This confirms the non-relativistic quark model Ansatz (i.e. = 2) as the 
preferred choice for this channel. The vector meson case is plotted in Fig. ^. Here, 
as expected, the signal is worse,! but the minimum is certainly not far oS n = 2. 
When these fits are applied to the nucleon, the data is more poorly behaved 
again (see Fig. However, it is clear that the minimum has shifted significantly 
to larger n values compared with the mesonic cases - it is now roughly between 
n — 5 and n = 10. 

We now study the behaviour of the Ansatz eq.(PD as a function of quark mass. 
In Fig. |10| the values for the pseudoscalar fits as a function of n are plotted 
for lighter quarks, corresponding to My f» 900MeV. (Recall that all other data 
in this work is for quark masses corresponding to My ^ l.lGeV - see Table 

^ In general, the pseudoscalar channel has the best signal to noise ratio. 
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Fig. 7. The value plotted against n for the pseudoscalar meson case for the fit 



corresponding to eq.(44)- The six lines are drawn as a guide for the eye only. They join 
points from the same simulations (see Table |^. 



[T|.) Comparing Fig. |T0| with Fig. |^ one can see a tendency for the minimum 
to decrease from n = 2 as the quark mass decreases. This fits nicely with our 
intuition - hghter quarks should eventually become relativistic which corresponds 
to n = (see Sec. 15.21). 



6 Discussion & Conclusion 



We have begun this work with a comprehensive study of the QCD Sum Rule 
Continuum Model as applied to lattice two-point correlation functions using the 
method first introduced by Leinweber |l|,0] . We have extended his work by includ- 
ing mesonic states, by fitting lattice data at several lattice spacings and by using 
two formulations of the lattice action. We have found that the QCD Sum Rule 
Continuum Model successfully fits the lattice data, and the quality of the fits is 
a significant improvement over the conventional "2-exponential" fits. However, 
while the results of the fits are successful, the values of the parameters in the fits 
are unphysical. This leads us to the conclusion (somewhat contrary to that of |jl[) 
that this model cannot self-consistently fit lattice data. 

We have then searched for a model that reproduces a similar functional form as 
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Fig. 8. As in Fig. but for the vector meson case. No points are plotted in cases where 
the fits did not converge. 




I 1 1 1 1 1 

2 4 6 8 10 



Fig. 9. As in Fig. but for the nucleon case. No points are plotted in cases where the 
fits did not converge. 
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Fig. 10. As in Fig. but for lighter quark masses corresponding to My k, 9QQMeV . No 
points are plotted in cases where the fits did not converge. 

the QCD Continuum Model, but which does not have the same restrictions on its 
parameters. Both the non-relativistic and relativistic quark models were studied. 
By using a fitting function which interpolated between the two, it was found that 
the non-relativistic quark model was qualitatively better that its relativistic coun- 
terpart for the quark mass studied. In the case of the pseudoscalar meson, these 
remarks are more quantitative. In addition, by studying lighter quark masses, the 
best fitting function was found to tend towards the relativistic form, as expected. 

Thus we have found a model which appears to correctly parametrise the role 
of the excited states. The functions that we have found can be widely used for 
Lattice Gauge Theories fits when the time separation in the correlation functions 
is forced to be small. In particular, it would be interesting to apply these functions 
to the case of the static case of the Heavy Quark Effective Theory on the lattice 
where ground state properties are poorly isolated (see, for example, [^]) and in 
glueball studies. 



Acknowledgement 

We wish to thank the APE collaboration for allowing us to use the lattice correla- 
tion function data presented here. We also thank Ian Drummond, Simon Hands, 



31 



Derek Leinweber, Vittorio Lubicz, Martin Liischer, Graham Shore and John Stack 
for useful discussions. CRA would like to acknowledge help from John Evans. 

S. C. thanks the Organising Committee for his participation at the 16th UK 
Institute at Swansea, where part of this work was carried out. 

This work was supported by the EC Contract "Computational Particle Physics" 
CHRX-CT92-0051, and by the EC Human Capital and Mobility Program, con- 
tracts ERBCHBICT941462 and ERBCHBGCT940665 and by a grant from the 
Nuffield Foundation. 



References 

[1] D.B.Leinweber, Phys. Rev. D51 (1995) 6369 
[2] D.B.Leinweber, Phys. Rev. D51 (1995) 6383 

[3] M.A. Shifman, A.I. Vainshtein and V.I. Zakharov, Nucl. Phys. B 147 (1979) 385, 
448 

[4] K.G.Wilson, Phys. Rev. DIO (1974) 2445, and in "New Phenomena in Subnuclear 
Physics", Plenum Press, New York (1977) (A.Zichichi ed.) 

[5] B.Sheikholeslami and R.Wohlert, Nucl.Phys. B259 (1985) 572. 

[6] M.-C. Chu, J. M. Grandy S. Huang and J. W. Negele, Phys. Rev. D48 (1993) 
3340. 

S.J.Hands, P.W.Stephenson and A.McKerrell (The UKQCD Collaboration), Phys. 
Rev. D51 (1995) 6394 

E. Shuryak, Rev. Mod. Phys. 65 (1993) 1. 



[9] C.R. Allton and S. Capitani, hep-lat /9709046 , to be published in Nucl. Phys. B 
(Proc. Suppl.). 

[10] APE Cohaboration, C.R. Allton et al, Phys. Lett. B 326 (1994) 295 

[11] APE Collaboration, C.R.Allton et al, Nucl. Phys. B (Proc. Suppl.) 42 (1995) 
385, and in preparation 

[12] M. Crisafulli et al, Phys. Lett. B 369 (1996) 325. 

[13] C.R. Allton, V. Gimenez, L. Giusti and F. Rapuano, Nucl. Phys. B489 (1997) 427. 
[14] UKQCD Collaboration, C.R.Allton et al, Phys. Rev. D 49 (1994) 474 
[15] B. L. loffe, Nucl. Phys. B188 (1981) 317. 



32 



[16] APE Collaboration (P. Bacilieri et al.), Nucl. Phys. B 317 (1989) 509. 

[17] JLQCD Collaboration (S. Aoki et al.), Nucl. Phys. B (Proc. Suppl.) 47 (1996) 
354. 

[18] Current Physics Studies and Comments, vol. 9, Quarkonia (ed. W. Buchmiiller) 
North Holland, 1992. 

[19] E. Eichten, Nucl. Phys. B (Proc. Suppl.) 20 (1991) 475, C. W. Bernard, Nucl. 
Phys. B (Proc. Suppl.) 34 (1994) 47, C. R. Allton, Nucl. Phys. B (Proc. 
Suppl.) 47 (1996) 31. 



33 



